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Abstract 

Wc study two l/Nc effects on the meson spectra by using the AdS/CFT correspondence 
where the l/N^ corrections from the chiral condensate and the quark density are controUed 
by the gravitational backreaction of the massive scalar field and U{1) gauge field respectively. 
The dual geometries with zero and nonzero current quark masses are obtained numerically. We 
discuss meson spectra and binding energy of heavy quarkonium with the subleading corrections 
in the hard wall model. 



*e-mail : cyong21@sogang.ac.kr 
^e-mail : bhl@sogang.ac.kr 
■'■e-mail : shin@theor.jinr.ru 



Contents 



1 Introduction 



1 



2 Asymptotic AdS background 



2 



3 Mesons with a zero niq 

3.1 Light meson spectra 

3.2 Binding energy of heavy quarkonium 



4 

4 
6 



4 Mesons with a nonzero niq 

4.1 Light meson spectra 

4.2 Binding energy of heavy quarkonium 



9 

9 
12 



5 Discussion 



12 



1 Introduction 

Holographic QCD provides descriptions of the strongly interacting regime of the gauge theory 
based on the anti-de Sitter/conformal field theory (AdS/CFT) correspondence [1, 2, 3]. The 
properties of the gauge theory including confinement, chiral symmetry breaking, supersymmetry 
reduction and glueballs are discussed in [4, 5, 6, 7, 8, 9, 10]. QCD-like models are proposed 
in [11, 12, 13, 14, 15, 16, 17, 18, 19, 20]. Baryons and chemical potential are considered in 
[11, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30]. 

One of the simplest bottom-up approach is the hard wall model [13, 14] where the confinement 
is realized by introducing an infrared cut-off in the AdS spacetime. The predictions of the hard 
wall model match measured values of meson masses within 10% error. 

AdS geometry is dual to the large Nc limit of U{Nc) gauge theory on the boundary, which is 
conformal. Pure AdS geometry can be considered as a UV fixed point of the holographic QCD. 
If we consider a subleading correction having a typical scale, the conformal symmetry should be 
broken. As a result, the dual geometry should be modified to include the gravitational backreaction 
of the bulk field dual to the subleading correction, which may improve the hard wall model. 

There are two important holographic subleading corrections. One is a correction from the 
hadronic medium and the other is a correction from the chiral condensate. The subleading correc- 
tion from the hadronic medium is the gravitational backreaction of local U{1) gauge fields. The 
Reissner-Nordstrom AdS black hole (RN AdS BH) and thermal charged AdS (tcAdS), which is 
the zero mass limit of the RN AdS BH, are proposed as the corresponding geometries for decon- 
finement and confinement phases respectively [28, 31, 32]. The Hawking-Page transition between 
the geometries, and the meson spectra and the binding energy of heavy quarkonium on the tcAdS 
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are discussed in [28, 31, 33, 34]. The subleading correction from the chiral condensate is the gravi- 
tational backreaction of a massive scalar field [35, 36, 37, 38]. In [35] the contribution of the scalar 
field to the Hawking-Page transition is investigated. In [36, 37], holographic models capturing the 
gravitational backreaction of the scalar field are constructed. In [38], the dual geometries with 
zero and nonzero quark masses are obtained by numerically solving the equations of motion of the 
bulk action including the scalar field. The light meson spectra and the binding energy of heavy 
quarkonium are discussed on both backgrounds. 

In this paper we study subleading corrections from the chiral condensate and the hadronic 
medium as an extension of [31, 33, 34, 38]. We numerically solve the equations of motion of the 
bulk action including a massive scalar field and U{1) gauge fields to obtain the dual geometries 
with zero and nonzero current quark masses. In QCD phenomenology, it is known that the chiral 
condensate gets reduced in the hadronic medium [39]. We constrain the dual geometries with a 
model-independent relation between the chiral condensate and the quark density. We calculate 
meson masses and binding energy of heavy quarkonium on the geometrical backgrounds. By 
comparing the results with [31, 33, 34, 38], we discuss the effects of the chiral condensate and the 
quark density. 

Organization of this paper is as follows. In section 2, the bulk action and the equations of 
motion of it for the dual geometries are discussed. In section 3 and 4, the equations of motion 
of the bulk action are solved with zero and nonzero current quark masses. The meson masses 
and the binding energy of heavy quarkonium are obtained on the asymptotic AdS geometries. In 
section 5, we summarize the results. 



2 Asymptotic AdS background 



The gravity action in the bulk is 

1 



d^xV-G 



2k2 
Tr 



(7^ - 2A) 



4^2 



piL) p(L) 



MN 



+ 



(2.1) 



where 



and A = We follow the convention that = and = [13]. We 

consider the case of Nf = 2 and Nc = 3. The superscripts (L) and (i?) are of SU (A^f )l x SU {Nf)ji 



flavor symmetry with F^j^'^^ 



of the complex scalar field is defined hy Dm^ = — ^^m''^ + 

1 



.{L,R) 



A 



{L,R) .{L,R) 



M ' 



. The covariant derivative 
We set 



(R) 

M 



(2.2) 



where the modulus of the complex scalar field is regarded as a background field giving the gravi- 
tational backreaction while 7r^{z) corresponding to the scalar meson is regarded as fluctuations. 
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The background geometries of our interest are obtained from a gravity action with the massive 
scalar field of the Lagrangian (2.1) and U{1) gauge interaction. The corresponding action is 



S = [ d^xV^l (7^ - 2A) - ^ r {dM^f + m^^^ 



(2,3) 



where (j){z) is the scalar field in (2.2). We choose an ansatz for the asymptotic AdS metric in the 
Feffer man- Graham coordinate as 

ds^ = {-F{z)dt^ + G{z)dx^ + dz^\ , (2.4) 

where R is the AdS radius. We take R = 1. For the AdS black hole without the chiral condensate 
and the quark density, the components of the metric (2.4) are given by 

(1 - Mzy 



F{z) 
G{z) 



l + Mz* 
1 + Mz*, 



(2.5) 



where M is the black hole mass, which is asymptotically the AdS space. In the AdS spacetime, 
the modulus of the scalar field (p(z) is 



(z) = uigZ + az^, 



(2.6) 



where ruq and a = (qq) are related to the current quark mass and the chiral condensate of QCD. 
The U{1) gauge fields in (2.3) are Aq = Ao{z) whereas Am = for M = 1, 2, 3, z. In the RN AdS 
BH or tcAdS, the time component vector field is 

Ao{z) = fi- Qz^ (2.7) 

where and Q are related to the chemical potential and the quark number density. 

Now, we generalize the metric (2.4) including the effects of the chiral condensate and the quark 
density. The Einstein equations and the equations of motion for (j){z) and Aq(z) are 

1 
4 
1 
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We solve the equations to obtain the geometries with zero and nonzero quark masses, and calculate 
the meson masses and the binding energy of heavy quarkonium on each geometry. 



3 Mesons with a zero rUq 
3.1 Light meson spectra 

We solve (2.8) with a zero quark mass, mq = 0, which corresponds to the chiral limit showing the 
spontaneous symmetry breaking effect caused by the chiral condensation. In the UV limit, the 
scalar field (p{z) and the vector field Aq{z) should be (2.6) and (2.7) respectively. We also use the 
fact that in the absence of the chiral condensate and the quark density, the solution should be 
(2.5). The perturbative solutions to (2.8) near the boundary, up to 0{z^'^) are 

F{z) = l-3Mz^ + — z^ + 4M^z^ + ^ -z^°, 



36^2 60^2 



48c/2 10 

^o(^) = ^ - Qz2 ^ ^Q^6 ^ - 3Qg'^'^% 8 _ M^Qz''. (3.1) 



We set M = 0, as we are interested in the meson spectra in the confining phase. From the above 
perturbative solutions, we can easily find the full numerical solutions depending on parameters, a 
and Q. Note that if we concentrate on the canonical ensemble is not important for investigating 
the meson spectra. 

The perturbation of the massive scalar field and the bulk gauge fields describe the pion, 
vector meson and axial-vector meson respectively. We transform perturbed gauge fields a^J'^^ of 
SU{Nf)L X SU{Nf)ji flavor symmetry to vector and axial-vector fields: 

= ^(ai/' + aif), Om = ^{afi^ - a[f). (3.2) 
The action describing the mesons is 



AS = / d^xV-G 



2 Ag^ 



(3.3) 



where the $ is as defined in (2.2) and fU^ and f^f}f are the field strengths of the gauge fields 
Vm and Qm- The Lagrangian is in the gauge where Vz = and = 0. The axial vector can be 
decomposed into a transverse component and a longitudinal component: 



a. 



«M + ^mX- (3.4) 
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We set Qq = as the Lorentz boost symmetry is not manifest. The equations of motion for Vi, ai, 
vr and x 

/111 \ 1 1 1 , 
d, i-F^G-^d.Vij +-F~2G2mlvi = 0, 

a, 0F^Gia,a,) - ^F-2G-2 (^V' - yn^l^ = 0, 

a, 0F-iGta,x) + 5''A'^i^-^Gi(7r - x) = 0, 

mldzX - 9^^4>^dzTT = 0, (3.5) 

where rriv, and m^r are the masses of vector meson, axial- vector meson and pion respectively. 
We solve (3.5) to obtain the meson masses in the background geometry numerically obtained from 
(3.1). 

The dual geometry corresponding to the ground state contains the chiral condensate as well 
as the quark density. The effects of the quark density and the chiral condensate on the meson 
spectra are discussed separately in [33, 38]. As the quark density increases, the meson mass and 
the oi-meson mass increase [33]. As the chiral condensate increases, the meson mass decreases 
whereas the ai-meson mass increases [38]. Thus it is worth investigating the effect of both XjN^ 
corrections to the meson spectra and which of them is more dominant. In phenomenology, the 
chiral condensate and the quark density are not independent. One of the known relations between 
them is 

— = 1-0.35—, (3.6) 

where pB the baryon number density and related to the quark number density with Q/Nc- 
The relation is considered to be model-independent, on the assumption that the pion-nucleon 
sigma term SttAt — 45MeV [39]. The quark density reduces the chiral condensate and it leads to 
the modification of the medium requiring non-trivial boundary conditions in the gravity set-up. 
Instead of finding such boundary conditions, we impose the relation (3.6) on the geometry (3.1). 

The quark density becomes a parameter. In our notation, Q/Qo = Pb/po with Qq = IGeV^ 
as a unit. We also choose 

do = (0.304GeV)^ z/^, = l/(0.3227GeV), (3.7) 

of which the physical aspects are discussed in [38]. We impose Dirichlet and Neumann boundary 
conditions at z = and z = zir respectively. The pion decay constant /tt is 

= _ ^M2) , (3.8, 

9 z 
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where a(0) is a solution of the second equation of (3.5) with = 0, satisfying the boundary 
conditions o(0) = 1 and dza{ziR) = [13]. The result is shown in Table 1, Figure 1 and Figure 2. 

The /9- meson mass increases as the quark density increases. This is qualitatively consistent 
with the meson spectra depending on the chiral condensate [38] and the quark density [33]. The 
ai-meson mass and pion mass, however, decrease with the quark density. This is consistent with 
the mass spectra depending on the chiral condensate [38], but contrary to the spectra depending 
on the quark density [33]. It shows that the effect of the chiral condensate dominates the effect 
of the quark density. 

The pion decay constant decreases with the quark density. From Gell-Mann-Oakes-Renner 
(GOR) relation 

f^ml = 2mga, (3.9) 
the deviation is A = f^m'^ since iriq = 0. The GOR relation is satisfied up to IQ-^^GeV''. 



Q(GeV^) 


mp(GeV) 


nia, (GeV) 


?n^(GeV) 


/.(GeV) 


A(GeV'*) 


0.00 


0.77581 


1.22166 


4.78271 X 10"*^ 


8.3073 X 10"^ 


1.5786 X 10"^^ 


0.01 


0.77583 


1.21970 


4.77344 X 10"*^ 


8.2946 X 10"^ 


1.5677 X 10^^^ 


0.02 


0.77591 


1.21779 


4.76435 X 10-*^ 


8.2818 X 10"^ 


1.5569 X 10"" 


0.03 


0.77603 


1.21594 


4.75546 X 10"*^ 


8.2691 X 10"^ 


1.5463 X 10-^^ 


0.04 


0.77621 


1.21415 


4.74675 X 10-*^ 


8.2565 X 10"^ 


1.5360 X 10~" 


0.05 


0.77643 


1.21242 


4.73825 X 10"*^ 


8.2439 X 10"^ 


1.5258 X 10"" 


0.06 


0.77671 


1.21074 


4.72994 X 10"*^ 


8.2313 X 10"^ 


1.5158 X 10"" 


0.07 


0.77703 


1.20913 


4.72182 X 10"*^ 


8.2188 X 10~^ 


1.5060 X 10"" 


0.08 


0.77740 


1.20757 


4.71390 X 10"*^ 


8.2063 X 10"^ 


1.4964 X 10"^^ 


0.09 


0.77782 


1.20607 


4.70618 X 10"*^ 


8.1939 X 10"^ 


1.4870 X 10^^^ 


0.10 


0.77829 


1.20463 


4.69867 X 10"*^ 


8.1815 X 10^^ 


1.4778 X 10^1^ 



Table 1: Masses of p-meson, ai-meson and pion. f^^ is the pion decay constant. A is the deviation 
from GOR relation, niq = 0, zjr = l/(0.3227GeV) and do = (0.304GeV)3. 

3.2 Binding energy of heavy quarkonium 

We study the binding energy of heavy quarkonium in the confining phase where the quark density 
is (3.6) and the asymptotic geometry is (3.1) with M = 0. We observe breaking of a fundamental 
string which connects two heavy quarks. It describes dissociation of a bound state of heavy quarks 
into two heavy-light quark bound states. The action for a fundamental string is 

MN, M.iU) 

ZTTa' J 
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mp[GeV] 




(b) (c) 
Figure 1: Meson masses with niq = 0. (a)p-meson (b)ai-meson (c)pion 




(a) (b) 
Figure 2: (a)pion decay constant (b)deviation from GOR relation, niq = 0. 



where a, b are the string worldsheet indices and Gmn is the metric (2.4). We take R = 1 and 
a' = l/27r. We choose a gauge condition and an ansatz for the coordinate z as 

T = t, ai = = X and z = z{x). (3-11) 
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The action becomes 

rT rr/2 1 

S= dt dx—JF(G + z'^). (3.12) 

Jo J-r/2 

Notice that we use the Fefferman-Graham coordinate and F and G are sohitions obtained by 
numerically solving (2.8), whose asymptotic expansions at the boundary are given in (3.1). The 
Hamiltonian density is 

-}{ = ^ (3 13) 

y/F{G + z'^y 

We define z = zq as the point where ff j^^^^^ = 0. The Hamiltonian density at z = zq becomes 

1 



-Ho = --VFhGo, (3.14) 

^0 

where Fq and Gq are the values of F and G at z = ^o- As Hamiltonian is conserved we can get 
the relation between r and zq from (3.13) and (3.14) as 

dz V » (3 ;^5) 







yJz^FG'^ - Z^FoGqG 

The kinetic energy of two heavy quarks can be ignored since we are considering a static configu- 
ration. The potential energy is obtained from (3.12) and (3.15) as 

V = 2 dz^ , , (3.16) 
Jo z^ V4^G^ - z^F^GoG 

The potential energy diverges since the heavy quarks have infinite masses at the boundary. We 
regularize the potential energy by subtracting the energy of two straight strings corresponding to 
two free heavy quarks. By choosing a gauge condition and an ansatz for the two strings as 

T = t, G\ = z and X = constant, (3-17) 

the energy of two static straight strings is 

rziR X r- 

Vo = 2 -^Vf. (3.18) 

Jo z 

The regularized binding energy of heavy quarkonium is 

E = 2 dz^ , ^ -2 (3.19) 

Jo V4PG^ - z^FoGoG Jo z^ 

The dissociation length increases with the quark density as shown in Table 2. It indicates that 
it takes more energy to produce a pair of heavy-light quark bound states as the quark density 
increases. By taking account of (3.6), this is consistent with the result of [38] that the dissociation 
length increases as the chiral condensate decreases. 
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g(GeV^) 


length with rrig = 0(GeV ^) 


0.00 


2.2217 


0.01 


2.2218 


0.02 


2.2220 


0.03 


2.2224 


0.04 


2.2228 


0.05 


2.2234 


0.06 


2.2242 


0.07 


2.2251 


0.08 


2.2261 


0.09 


2.2272 


0.10 


2.2285 



Table 2: Dissociation length, rrig = 0, zm = l/(0.3227GeV) and ao = (0.304GeV)^. 



4 Mesons with a nonzero ijiq 
4.1 Light meson spectra 

We solve (2.8) with a nonzero quark mass, nig = 0.002383GeV 
symmetry explicitly. The asymptotic solutions are 



which breaks the chiral 



F{z) 



G{z) 



1 

12 
1 



17280 



1 

144 



K^m^z^ + TTrii^^fTT'g ~ ISK^rUga — 3^^771^ log 



— + 65K^m^ - 396«;'^m^cj - 2592^^^ 



31104 \ r 

- 66K^m^ log z - 864K''m^cj log z - 72K''m^(log z)'^ ) z*" + 



1 

12 

1 



1 



-t -*-222, ^ / d d -in2 o44i \4 

1 — —K m^z H — ttI'^ '^g ~ lo/^ mgcr — 3k log Z)z 



144 



3456- 



^ + 65K'^m^ - 396K^m^f7 - 2592^^^ 



31104 \ g 

- mK^ml log z - ^UK^mla log z - 72K^m^(log zf ] z^ + 



niaZ + fcr + -K^ml logz]z^ + — (-13K''m^ + lA4.K^mla + 24K''m^ \ogz)t^ 

y ^ g g o y 576 go/ 

1 / 2"^g'3^ 



+ 



62208 



1728k^ 



9' 



+ llAK^'m'g log z + 3456fi;* m^o- log z + 288k'' (log zY]z' + 

9 




+ 'iUQn^ml log z + lllSQti^mla log z + UAQK^ml{\og zf\z^ + '"- 



(4.1) 



To investigate the meson masses, we solve (3.5) on the full numerical solutions obtained from 
(4.1). We also use the relation (3.6) with Qq as a unit and the values of the chiral condensate and 
IR cutoff (3.7). The result is shown in Table 3, Figure 3 and Figure 4. 



Q(GeV^) 


mp(GeV) 




m^(GeV) 


/.(GeV) 


A(GeV^) 


0.00 


0.77580 


1.23056 


0.139617 


8.4645 X 10"^ 


5.7645 X 10"*^ 


0.01 


0.77583 


1.22860 


0.139590 


8.4519 X 10^^ 


5.7631 X 10"*^ 


0.02 


0.77590 


1.22669 


0.139568 


8.4394 X 10-^ 


5.7756 X 10"« 


0.03 


0.77603 


1.22484 


0.139551 


8.4269 X 10^^ 


5.8017 X 10"*^ 


0.04 


0.77620 


1.22305 


0.139541 


8.4145 X 10^^ 


5.8416 X 10"*^ 


0.05 


0.77643 


1.22132 


0.139536 


8.4021 X 10"^ 


5.8956 X 10"*^ 


0.06 


0.77670 


1.21965 


0.139538 


8.3897 X 10-^ 


5.9634 X 10-*^ 


0.07 


0.77702 


1.21803 


0.139545 


8.3774 X 10^^ 


6.0454 X 10"*^ 


0.08 


0.77740 


1.21648 


0.139559 


8.3652 X 10^^ 


6.1418 X 10"*^ 


0.09 


0.77782 


1.21498 


0.139579 


8.3530 X 10^^ 


6.2523 X 10"*^ 


0.10 


0.77829 


1.21354 


0.139606 


8.3408 X 10~^ 


6.3772 X 10-*^ 



Table 3: Masses of /9-meson, ai-meson and pion. /^r is the pion decay constant. A is the deviation 
from GOR relation, niq = 0.002383GeV, zm = l/(0.3227GeV) and ctq = (0.304GeV)3. 

The /9-meson mass increases whereas the ai-meson mass decreases as the quark density in- 
creases. The patterns are the same as the results with a zero quark mass as observed in Section 
3.1. The ai-meson mass decreases with the quark density in both cases of zero and nonzero quark 
masses. By taking the relation (3.6) into account, the spectrum is consistent with the result of 
[38] that the ai-meson mass increases with the chiral condensate, but contrary to the result of 
[33] that the ai-meson mass increases with the quark density. It shows that the effect of the chiral 
condensate on the mass spectra dominates the effect of the quark density. A distinguishing feature 
is that the spectrum of pion exhibits a critical point in the case of a nonzero quark mass whereas 
it decreases monotonically with the quark density in the case of a zero quark mass as observed 
in Section 3.1. At the low quark density, the chiral condensate is high so that the quark mass 
is negligible. Thus the pion mass spectrum is similar to the case of a zero quark mass. At the 
high quark density, however, the chiral condensate is small. The effect of the quark mass becomes 
comparable with the chiral condensate near the critical point. Above the critical point, the effect 
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(b) (c) 
Figure 3: Meson masses with ruq = 0.002383GeV. (a)p-meson (b)ai-meson (c)pion 



/jGeV] 




g[GeV 

0,10 




0,10 



Figure 4: (a)pion decay constant (b)deviation from GOR relation. rUq = 0.002383GeV. 



of the quark mass, which is constant, dominates the chiral condensate so that the quark density 
mainly contributes to the medium. The pion mass therefore increases with the quark density as 
it is observed in [33]. 

The pion decay constant decreases as the quark density increases. The deviation from the 
GOR relation (3.9) is A = f^ml - 2mqa. The GOR relation is satisfied up to IQ-^GeV^. 
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4.2 Binding energy of heavy quarkonium 

We study the binding energy of heavy quarkonium in the confining phase. We numerically solve 
the equations (2.8) with (4.1) as asymptotic expansions of the metric components at the boundary. 
The binding energy of heavy quarkonium is (3.19) and the dissociation occurs when E = 2mq. 
The result is shown in Table 4. 

The quark mass is small so that there is no big difference between the cases of zero and nonzero 
quark masses as shown in Table 2 and Table 4. The dissociation length with a nonzero quark 
mass is slightly longer. It implies that the quark mass shifts the binding energy a little bit higher. 



g(GeV^) 


length with rUg / 0(GeV"^) 


0.00 


2.2382 


0.01 


2.2383 


0.02 


2.2385 


0.03 


2.2389 


0.04 


2.2394 


0.05 


2.2400 


0.06 


2.2407 


0.07 


2.2416 


0.08 


2.2426 


0.09 


2.2438 


0.10 


2.2450 



Table 4: Dissociation length. rUg = 0.002383GeV, zjr = l/(0.3227GeV) and ctq = (0.304GeV)3. 



5 Discussion 

We have studied sub leading 1/Nc corrections from the chiral condensate and the quark density 
to the meson spectra and the binding energy of heavy quarkonium. We have considered the 
gravitational backreaction of a massive scalar field, which corresponds to the current quark mass 
and the chiral condensate, and the time component of U{1) gauge fields, which corresponds to the 
chemical potential and the quark density. The geometries are numerically solved, separately with 
zero and nonzero current quark masses. The geometries are constrained by a model-independent 
relation in QCD phenomenology that the chiral condensate gets reduced linearly with the quark 
density. The meson masses and binding energy of heavy quarkonium are calculated on each 
geometry. 

With a zero quark mass, the p-meson mass increases as the quark density increases. This is 
consistent with the results of [33, 38], where the meson spectra are studied with the quark density 
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and the chiral condensate as parameters separately. The ai-meson mass and pion mass decrease 
with the quark density. The results are consistent with the spectra with the chiral condensate, 
but contrary to the spectra with the quark density. 

With a nonzero quark mass, the spectra of the p-meson and ai-meson present the same pattern 
as the spectra with a zero quark mass. The oi-meson mass decreases with the quark density in 
both cases of zero and nonzero quark masses. By taking account of the fact that the chiral 
condensate gets reduced with the quark density, this is consistent with the result of [38] that the 
ai-meson mass increases with the chiral condensate, but contrary to the result of [33] that ai- 
meson mass increases with the quark density. It indicates that the effect of the chiral condensate 
on the mass spectra dominates the effect of the quark density. The mass spectrum of pion exhibits 
a critical point. Below the critical point, where the quark density is low, the chiral condensate is 
high so the quark mass is negligible. Thus the pion mass decreases with the quark density as it 
does with a zero quark mass. Above the critical point, where the quark density is high and the 
chiral condensate is low, the effect of the constant quark mass dominates the chiral condensate. 
As the quark density becomes the main contribution to the medium, the pion mass increases with 
the quark density as observed in [33] . 

The dissociation length of heavy quarkonium increases with the quark density for both cases 
of zero and nonzero quark masses. It shows that it requires more energy to produce a pair of 
heavy-light quark bound states as the quark density increases. This result is consistent with the 
dissociation length depending on the chiral condensate [38]. The dissociation length is slightly 
longer with a nonzero quark mass. It indicates that the quark mass shifts the binding energy 
higher. 
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